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The plasmonic resonances of nanostructured silver films produce exceptional surface enhancement, 
enabling reproducible single-molecule Raman scattering measurements. Supporting a broad range of 
plasmonic resonances, these disordered systems are difficult to investigate with conventional far-field 
spectroscopy. Here, we use nonlinear excitation spectroscopy and polarization anisotropy of single optical 
hot spots of supercontinuum generation to track the transformation of these plasmon modes as the 
mesoscopic structure is tuned from a film of discrete nanoparticles to a semicontinuous layer of aggregated 
particles. We demonstrate how hot spot formation from diffractively-coupled nanoparticles with broad 
spectral resonances transitions to that from spatially delocalized surface plasmon excitations, exhibiting 
multiple excitation resonances as narrow as 13 meV. Photon-localization microscopy reveals that the 
delocalized plasmons are capable of focusing multiple narrow radiation bands over a broadband range to the 
same spatial region within 6 nm, underscoring the existence of novel plasmonic nanoresonators embedded 
in highly disordered systems. 

Localized surface plasmon resonances (LSPRs), the resonant oscillations of conduction electrons of nano- 
structured metallic systems, enable access to experiments involving macroscopic electronic coherence 1 ' 2 . 
Since this resonance is often found at visible and infrared frequencies, plasmonic nanostructures provide an 
attractive optoelectronic platform that enables the near-field control of light on dimensions that fall far below the 
diffraction limit 3 " 5 . For example, metallic nanostructures are capable of focusing incident electromagnetic radi- 
ation to minute, nanoscale regions, termed "hot spots," where the electric field can be enhanced by several orders 
of magnitude over that of the incident field. Within the hot spots, the interaction between light and matter can be 
significantly amplified giving rise to surface-enhanced optical phenomena 6 . One of the foremost exploitations of 
this effect is surface-enhanced Raman scattering (SERS) where the interactions between light and the vibrational 
modes of matter within the hot spot can be magnified, enabling Raman scattering measurements of a single 
molecule. Such a feat necessitates overall enhancement factors approaching 12 orders of magnitude 7 " 10 . Beyond 
SERS, plasmonic systems have also been employed in a wide range of nanophotonic applications such as near- 
field imaging 11 and transmission spectroscopy 12 using highly localized nonlinear light sources; nanoscale lasers 13 ; 
and ultra-sensitive analysis of dielectric environments 14 . 

Significant effort is placed in designing metallic nanostructures, and because of their electrodynamic nature, 
the parameter space that is available for engineering their plasmonic resonances is expansive and permits a wide 
range of possibilities 1415 . For a single metallic nanoparticle, shape, size and composition are the primary toolset, 
and by utilizing the near-field and far-field coupling 516 , the parameter space is greatly expanded for multi-particle 
arrangements. In contrast to carefully designed structures, one of the most attractive and fascinating multi- 
particle systems consists of random metal -dielectric films 17 " 20 . The naturally arising structural intricacy of these 
films yields plasmonic resonances that result from a complex hybridization between spatially localized modes and 
delocalized modes that extend over mesoscopic dimensions of the film 141719 2126 . Remarkably large enhancement 
factors are found within the resulting hot spots, and it is with these structures that numerous single molecule SERS 
studies have been performed 910,27 . Further illustrating the advantageous nature of these plasmon resonances, 
evolutionary computational algorithms that are designed to optimize electromagnetic enhancement of a metal 
nanostructure have been shown to generate a disordered structure as opposed to an ideal dipolar antenna 18 . A 
thorough understanding of electronic coherence and plasmon resonances in these disordered systems has 
potential to provide novel approaches to the design of highly optimized nanostructures and motivates theoretical 
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and experimental studies of the electronic structure-property rela- 
tionships that govern the behaviour of these modes. 

In this work, we use the nonlinearly excited supercontinuum emis- 
sion to reveal the plasmonic modes of nanostructured silver films 
produced with the Tollens reaction 27 as the morphology is varied 
from discrete nanoparticles to a rough, semicontinuous metal layer, 
using different growth times 28 . The nearly background-free emission 
originates from discrete diffraction -limited regions, and its exact 
nature and origin remain poorly understood. However, since the 
emission results from plasmonic hot spots 20,28 , the excitation process 
provides a facile, noninvasive probe to study the intricate plasmon 
resonances of the disordered silver nanoparticle system. The spectral 
characteristics of the plasmon modes are determined from nonlinear 
excitation spectroscopy of single hot spots and exhibit a remarkable 
change from single-peaked, broad resonances in the low-coverage 
films to narrow, multi-peaked resonances in high-coverage films. 
Through measurements of excitation polarization anisotropy of sin- 
gle hot spots, we can discern signatures of delocalized modes result- 
ing from the diffractive (dipolar) coupling in the low-coverage films 
as opposed to the (non- diffractive) derealization of the plasmon 
modes in the extended metal structures of the high-coverage films. 
The surface coverage at which significant plasmon derealization is 
expected is marked by a drop in the number of nonlinear hot spots 
that are excitable with an electric field oriented parallel to the sample 
surface, suggesting a conceptual connection to propagating surface 
plasmons in a smooth metal film 29 . At the highest surface coverage 
where the modes are expected to be maximally delocalized, the films 



consist of multi-resonant hot spots that can only be excited with 
electric fields that are normal to the sample surface. Using photon- 
localization microscopy 30 , we provide evidence that these multi- 
peaked resonances in the semicontinuous films belong to one and 
the same hot spot. 

Results 

Polarization anisotropy and excitation spectroscopy of single hot 
spots. The nonlinear optical excitation measurements are illustrated 
in Fig. 1 for a high- coverage, semicontinuous silver film (growth time 
255 s). Figure 1(a) shows the illumination scheme used for the 
nonlinear optical microscopy where linearly polarized, pulsed, 
infrared (IR) radiation illuminated the sample at an oblique angle 
of incidence. The polarization of the nonlinear excitation is described 
by q> ki the angle that is formed between the horizontal direction and 
the electric field vector of the optical excitation. Under p-polarization 
(<Pk = 0°, 180°), the electric field lies in the plane of incidence and 
consists of one component that is parallel to the sample plane and an 
orthogonal component that is normal to the sample plane. For s- 
polarization (q? k = 90°, 270°) the electric field is entirely parallel to 
the sample plane (see the Supplementary Information for a complete 
discussion of the polarization). Inset in Fig. 1(a) is a representative 
SEM micrograph of the semicontinuous film showing that the 
structure is marked by a dense distribution of nanoscale voids and 
protrusions. Alongside the SEM image, an optical micrograph of the 
film under excitation at 1.78 eV shows that the upconverted 
supercontinuum arises from discrete diffraction -limited spots, 
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Figure 1 | Excitation spectroscopy and polarization anisotropy of single nonlinear hot spots in a high- coverage, semicontinuous silver film grown 
following the Tollens mirror reaction, (a) Under illumination of obliquely incident, pulsed infrared radiation, the film exhibits surface-enhanced 
upconverted supercontinuum emission from discrete, diffraction-limited spots. Typical SEM and fluorescence micrographs are shown, (b) Excitation 
spectrum of a single hot spot where the emission intensity is plotted as a function of the IR photon energy. The black line serves as a guide to the eye. The 
hot spot exhibits multiple resonances that are each characterized by their FWHM, AE. (c) Excitation polarization anisotropy measurements on single hot 
spots. The emission intensity is recorded as a function of the angle of the electric field vector from the horizontal, (p^ defined in the cartoon, and can be 
fitted with a cos 4 function (black line). The coupling to polarized excitation is characterized by the polarization anisotropy value M, defined in the inset, 
and the in-plane dipole angle, <^dipoie> which is determined from the angle of maximum emission (p™*. 
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which we refer to as "nonlinear hot spots." Temporal and spectral 
analysis of the emission from these nonlinear hot spots 28 reveals that 
it is quasi-instantaneous with a lifetime of less than 5 ps and is 
comprised of a broadband continuum with a small amount of 
fluctuating second-harmonic radiation 31 , which can be neglected in 
spectrally integrated measurements. Figure 1(b) shows an example 
excitation spectrum of a single nonlinear hot spot from a semi- 
continuous film where the emission intensity is recorded as a 
function of excitation energy. Each excitation resonance was 
characterized by its full width at half maximum (FWHM), AE. 
Additionally, we preformed excitation polarization anisotropy 
measurements on single hot spots. A sample dataset from such a 
measurement is shown in Fig. 1(c), where the intensity of the 
continuum emission is plotted as a function of the excitation 
polarization in terms of (p^. From such traces, the polarization- 
dependent coupling between the hot spot and IR excitation is 
summarized in two extracted values: the angle of the excitation 
dipole in the sample plane, dipole ' which can be calculated from 
the corresponding angle of maximum emission intensity, <^™ ax , and 
the anisotropy value, M, which quantifies the emission modulation 
depth 32 . 

Figure 2 shows the distributions of the resonance widths, AE, 
(Column A), excitation anisotropy values, M, (Column B) and in- 
plane dipole angles, ^ dipo i e , (Column C) for three different films 
(Rows I— III) with distinct surface morphologies as summarized in 
the SEM micrographs inset in Column C. Progressing downward 
from the topmost row, the surface coverage of the three films 
increases from well- separated nanoparticles deposited with a short 
growth time (45 s; Row I) over interconnected silver islands (inter- 
mediate growth time of 135 s; Row II) to a semicontinuous structure 
with a dense distribution of nanoparticles (long growth time of 255 s; 
Row III). In previous studies, it was shown that the growth time 
provides a reliable parameter for reproducible structures. The films 
of interconnected islands were found to yield maximal linear and 



nonlinear emission 28 as well as optimal SERS performance 33 . All 
three distributions of the excitation properties display a clear 
dependence on film morphology. For each film, the distribution of 
resonance widths is well described by a single Gaussian centred 
around an average value, AE. The nonlinear hot spots in the low- 
coverage film (Row I) exhibit mostly single-peaked resonances with a 
broad distribution of widths centred on AE = 145 meV with a 
standard deviation (SD) of 50 meV. An increase of growth time 
(Row II) leads to a narrower distribution of resonance widths centred 
on a notably smaller average value, AE = 105 meV with a SD of 
33 meV. This narrowing is more dramatic in the semicontinuous 
films (Row III) which show a distribution of widths centred on a 
significantly smaller average value, AE = 39 meV with a SD of 
13 meV. This pronounced decrease in resonance width as the surface 
coverage increases is accompanied by a transformation from single- 
peaked excitation resonances to multi-peaked resonances as seen in 
representative excitation spectra inset. 

The film morphology also affects how the nonlinear hot spots 
selectively couple to the laser polarization. As seen in Column B, 
the low-coverage film (Row I) yields a bimodal distribution of aniso- 
tropy values with a strongly-polarized population centred on M = 
0.95 and a less-polarized population around M = 0.65. The overall 
distribution is well -described by two super-imposed Gaussians. The 
relative populations of the strongly-polarized species to the less- 
polarized species are characterized by the heights of the respective 
Gaussian peaks in the histograms, P M >o.9 and P M <o.9- As the surface 
coverage is increased, excitation of the nonlinear hot spots becomes 
increasingly more anisotropic in terms of polarization, and the popu- 
lation of less-polarized hot spots, P M <o.9> vanishes to nearly unde- 
tectable levels in both the film of interconnected islands (Row II) and 
the semicontinuous film (Row III). 

A similar trend is also observed in the distributions of the in -plane 
dipole angles (Column C). The low-coverage silver film (Row I) yields 
a histogram that consists of a pronounced peak at 0° superimposed on 
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Figure 2 | Nonlinear excitation properties of single hot spots of Tollens silver films in three distinct regimes of surface coverage and morphology. For 

each film, the distribution of resonance widths (Column A), polarization anisotropy values (Column B) and in-plane dipole angles (Column C) are 
shown. Representative excitation spectra of single hot spots are inset in Column A; SEM micrographs of the substrates are given in Column C. At low 
surface coverage (Row I) the films consist of well- separated silver nanoparticles. Increased surface coverage leads to films that are made up of 
interconnected silver islands (Row II) and are ultimately semicontinuous (Row III). The distribution P of anisotropy values is bimodal, corresponding to 
unpolarized and polarized hot spot excitation. The distribution of phase angles arises from hot spots which can be excited by either p- or s-polarized light, 
the number of which is described by N. 
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top of a uniform distribution that spans all angles resulting in a 
positive baseline offset of the peak. This dual- featured scatter of angles 
in itself is surprising. For an isotropic arrangement of dipoles in two 
or three dimensions, the strong peak at 0° is unexpected (see 
Supplementary Information), and leads us to conclude that it must 
arise due to a preferential response of the nonlinear hot spots to p- 
polarized excitation where the electric field component of the incident 
radiation that is normal to the sample plane is maximized. The uni- 
form distribution that spans all angles, on the other hand, is expected 
for a random arrangement of dipoles in a two-dimensional plane (see 
Supplementary Information) and is accordingly attributed to hot spots 
that are maximally excited with in-plane electric fields. Thus, the low- 
coverage film consists of two types of nonlinear hot spots: p-polarized 
hot spots that respond maximally to an electric field that is normal to 
the sample plane; and s-polarized hot spots that respond maximally to 
an electric field that is parallel to the plane. Analogous to the aniso- 
tropy (M) distributions, these two populations are characterized by 
the height of the baseline offset for the s-polarized population (N s _ po i) 
and the height of the Gaussian peak for the p-polarized population 
(Np-poi) in the histograms in Column C. The values of M and ^ Qax for 
a single hot spot were found to be uncorrelated indicating that these 
two populations are independent of one another. So, in addition to the 
pronounced narrowing of resonance widths, increased surface cov- 
erage of the silver films also reduces the number of less-polarized and 
s-polarized hot spots leaving only highly-polarized hot spots that 
maximally couple to p-polarized excitation. 

Using the average resonance width (AE), the ratio between the 
less-polarized and strongly-polarized hot spots (Pm<o.9^m>o.9) an d 
the ratio between the populations of s-polarized andp-polarized hot 
spots (N s _ po i/Np_ po i), we summarize the progression of the excitation 
characteristics with increasing surface coverage in Fig. 3. In order to 
provide a more quantitative description of film morphology than 
growth time, we approximated surface coverage by computing the 
two-dimensional filling factor from SEM micrographs for eight dif- 
ferent growth times from 45 s to 255 s as shown in Fig. 3(a). The 
filling factor increases almost linearly from ~0 to a saturation value 
close to unity. The progression of the hot spot excitation properties is 
quite pronounced. As shown in Fig. 3(b), the average resonance 
width narrows monotonically with increased growth time. On the 
other hand, the relative populations of the less-polarized hot spots to 
strongly-polarized hot spots, as well as the ratio of the populations of 
s-polarized hot spots to p-polarized hot spots, show a pronounced 
step-like dependence on growth time as seen in Figs. 3(c) and 3(d), 
respectively. The populations of less-polarized (P M <o.9) an d s ~ 
polarized (N s _ po i) hot spots drop to nearly undetectable levels upon 
increasing the surface coverage of the film from a filling factor of 0.4 
(growth time: 105 s) to 0.6 (growth time: 135 s). This range of sur- 
face coverage is close to the percolation threshold of thermally depos- 
ited silver films 21 ' 24,25 . 

From Fig. 3, the nature of the plasmonic resonances is clearly 
distinct in the different regimes of surface coverage. When the film 
consists of well-separated silver islands, the resonances are domi- 
nated by the plasmon modes of interacting particles as is particularly 
evident in the broad resonance widths. In this regime, it is expected 
that the most pronounced hot spots are formed in the gaps of closely- 
spaced collections of particles, as in the case of simple arrangements 
such as nanoparticle dimers and trimers 34 . Such structures are 
expected to be maximally excited by electric fields that lie parallel 
to the sample surface 35 . Indeed, such arrangements are readily 
observed in our polarization data as the population of s-polarized 
hot spots. However, a larger population of p-polarized hot spots 
exists in these films and indicates that the majority of hot spots are 
maximally excited by an electric field that is normal to the sample 
surface. We attribute this pronounced preference in electric field 
vector orientation to diffractive coupling between the nanoparticles 
since light that is polarized normal to the sample plane maximizes the 
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Figure 3 | Transformation of plasmon modes with surface coverage as 
probed by nonlinear excitation of single hot spots, (a) The two- 
dimensional filling factor is used to describe surface coverage, showing an 
approximately monotonic increase with growth time. For low-coverage 
films, the plasmon modes are dictated by LSPRs of silver nanoparticles that 
are diffractively coupled to one- another, while the higher- coverage films 
are characterized by delocalized modes as is illustrated in the insets, (b) The 
average resonance width plotted on top of the corresponding distribution 
of data points shows a clear narrowing with growth time. The small data 
points describe the underlying distribution of resonance widths at the 
specified growth time and the horizontal breadth of the points reflects the 
occurrence frequency. The ratio between unpolarized to polarized (c) and 
s-polarized to p-polarized (d) hot spots shows a step -like drop with 
increasing growth time between 105 s and 135 s. 

out-of-plane polarization of the individual particles and concomi- 
tantly the in -plane diffractive coupling in the system 36 ' 37 . The result- 
ing picture for the low- coverage films is an intricate electromagnetic 
confinement scheme that combines both the localized nature of 
the near-field coupled LSPR modes of closely- spaced nanoparticles 
with the delocalized nature of the longer-range diffractive coupling, 
giving rise to electromagnetic localization reminiscent of disordered 
media 38 . 

As the silver coverage increases and larger, interconnected islands 
and ultimately semicontinuous films are formed, the plasmonic 
excitations transform from the LSPR modes of the individual 
nanoparticles to hybridized localized and delocalized modes, a well 
established phenomenon in semicontinuous metal films 1 ' 4,21 " 26 . 
Signatures of this transformation in the macroscopic optical prop- 
erties of such films can be treated using an effective medium theory 
such as the Maxwell-Garnett approximation 39 . However, the spatial 
averaging of the local field distribution in such theoretical treatments 
limits their applicability to high resolution analysis of optical hot spot 
formation 40 . Past measurements on thermally- deposited semi- 
continuous metal films using near-field microscopy 24 ' 25 , far-field 
spectroscopy 21 and fluorescence lifetime microscopy 26 have corre- 
lated the percolation threshold with the maximum level of hybrid- 
ization between localized and delocalized plasmonic modes. From 
the step-like behaviour near the percolation threshold that is seen in 
Fig. 3, we propose that the polarization anisotropy serves as an 
almost binary indicator of delocalized modes in these films: a 
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plasmonic mode with traits of derealization is maximally excited 
with p -polarized radiation. Such a condition on the excitation polar- 
ization is analogous to quintessential delocalized plasmon modes: 
propagating surface plasmons in smooth metal films, which neces- 
sitate a non-zero electric field that is normal to the metal surface 29 . 

Increased derealization is accompanied by a transformation of 
the excitation spectra. The observed resonances in the films of high- 
est surface coverage (i.e., the excitation spectrum in Fig. 2, Column A, 
Row III) consist of multiple narrow peaks, differing notably to the 
broad, single-peaked resonances that are observed in the low- 
coverage films (cf. Row I). In contrast to the step-like behaviour of 
the polarization data, the resonance widths narrow continuously 
with increased surface coverage, conceptually bridging the broad 
resonant scattering spectra of individual metal nanoparticles 41 ' 42 to 
narrow resonant near-field scattering spectra of large colloidal aggre- 
gates 17 . At the single hot spot level, devoid of ensemble broadening, 
we have observed resonances as narrow as 13 meV (255 s growth 
time). As spectral width and mode lifetime are inversely propor- 
tional, the narrow resonances suggest that the dephasing time of 
these modes can be over an order of magnitude larger than what is 
observed in single silver nanoparticles 41 ' 42 , qualitatively confirming 
an analogous conclusion derived from near-field scattering spectro- 
scopy 17 . We propose that spectral narrowing and the associated 
increase in lifetime is yet another signature of plasmonic derealiza- 
tion, as is particularly evident with the progression of these features 
with increasing surface coverage. In addition to this lifetime increase, 
the narrow resonances may have substantial practical appeal in ana- 
lytical applications where sensitivity is inversely proportional to spec- 
tral width 43 . Indeed, extraordinarily narrow resonances have been 
observed in well- engineered systems such as periodic arrays of nano- 
particles 44 " 46 . Being significantly more disordered, however, the 
rough metal films presented here are clearly distinct, and the exist- 
ence of sharp resonances in these systems has previously only been 
explored in the context of linear Rayleigh scattering in near-field 
microscopy 1719 and photoelectron emission microscopy (PEEM) 47 
of single hot spots, which do not necessarily correlate with the under- 
lying plasmonic excitations. In contrast, we report here directly on 
narrow plasmonic resonances detected in the optical far-field. 

At the highest coverage, single hot spots exhibit multiple reso- 
nances, suggesting that the film contains plasmonic structures capable 
of efficiently focusing several discrete bands of electromagnetic energy 
to the same spatial region (i.e., the same hot spot). However, a dense 
distribution of nanoscale features exists (cf. inset in Fig. 2, Column C, 
Row II) on scales well below the diffraction limit of our microscope 
(—500 nm), and the anticipated size of a hot spot is —15 nm 48 . Thus, 
extreme care must be taken when assigning the distinct resonances to 
the same spatial region. Even though near- field Rayleigh scattering of 
colloidal silver clusters has shown that a 100 nm region can contain 
multiple spectral resonances 1719 , this length scale is nearly an order of 
magnitude larger than the hot spot. Since the spatial distribution of 
near-field intensity depends dramatically on incident photon 
energy 41719 , the existence of hot spots with multiple discrete electro- 
magnetic resonances remains to be unambiguously established, neces- 
sitating spatial resolutions of order 10 nm. 

Localization microscopy of single multi-resonant nonlinear hot 
spots. To noninvasively probe the spatial extent of the multiple 
resonances of single hot spots in Figure 4, we employed photon 
localization 30 , a super- resolution microscopy technique that has 
been used in plasmonics to map the intensity distribution in the 
near-field regime around bow-tie antennae 49 , SERS hot spots in 
silver nanoparticle aggregates 50,51 and hot spots that enhance the 
fluorescence of single molecules on rough silver films 48 . As shown 
in Figure 4(a), the emission centroid of the hot spot is extracted by 
fitting a two-dimensional Gaussian to the diffraction-limited spot. 
Using these centroid coordinates, a single emitter can be localized 
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Figure 4 | Super-resolution localization microscopy of single multi- 
resonant nonlinear hot spots, (a) The centroid of the diffraction-limited 
hot spot image is extracted from a two-dimensional Gaussian fit. 
(b) Excitation spectrum of a single hot spot with corresponding 
microscope images. The hot spot centroid at each resonance was recorded 
20 times in an alternating fashion, (c) Resulting spatial localization of the 
hot spot at each resonance with respect to the average value, illustrating 
that the centroid at each resonance is localized within a region of ~6 nm. 
Repeating this measurement for 5 additional hot spots that included 25 
distinct resonances provides a distribution of spatial shifts of the centroid 
position between resonances and shows that multi-resonant hot spots are 
localized within an average region of ~ 10 nm. 

with a resolution limited only by the noise in the measurement, 
readily reaching nanometer precision 30,48 . The excitation spectrum 
and the corresponding microscope images of the hot spot at each 
resonance energy are shown in Fig. 4(b). The diffraction-limited 
emission profile of the hot spot is well described by a two- 
dimensional Gaussian function at each resonance. Furthermore, 
even though the emission spectrum varies from spot to spot 28 , it is 
virtually identical at each resonance of a single hot spot and over 10X 
broader (>400 meV) than the individual excitation resonances (see 
the Supplementary Information). Figure 4(c) shows the relative 
locations of the centroid at each resonance. Thus, for this hot spot, 
the centroid at all five resonances is localized to the same region 
within ~6 nm. Analogous measurements were repeated for 5 
separate hot spots comprising 25 distinct spectral resonances, and 
the corresponding distribution of spatial shifts, Ar, between 
resonances within a hot spot is shown in Fig. 4(d). All of the 
recorded shifts fall below 30 nm with a large number of shifts 
below 15 nm. We stress that these shifts span a range of 
—500 meV in excitation energy (over 200 nm in wavelength) and 
that such highly-localized emission is an unlikely scenario if the 
emitter were spatially extended. While the precise mechanism of 
nonlinear supercontinuum generation remains elusive, the high 
level of localization suggests the presence of a localized polarizable 
species within the hot spot 28 ' 31 ' 50 . A silver or silver oxide nanocluster is 
an obvious candidate for such a species 52 . By using this localized 
emitter as a probe of both localized and delocalized plasmon 
modes that are responsible for the nonlinear excitation process, we 
have identified evidence of naturally-grown plasmonic structures 
capable of producing multi- resonant broadband hot spots. 



Discussion 

We have shown that the excitation properties of surface-enhanced 
nonlinear supercontinuum hot spots in rough silver films grown by 
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the Tollens reaction reflect the signatures of the plasmonic reso- 
nances within these systems. For the low- coverage films that consist 
of well-separated silver islands, excitation spectroscopy and polar- 
ization anisotropy reveal that the hot spots form by diffractive coup- 
ling of the LSPRs of the individual nanoparticles. Signatures of 
non-diffractive delocalized plasmon resonances in the higher- 
coverage films that are near the percolation threshold are found in 
spectral narrowing of the resonance and a reduction in number of 
both hot spots that can be excited with electric fields oriented parallel 
to the sample plane, and hot spots that have reduced excitation 
polarization anisotropy. We conclude that delocalized modes have 
narrow resonances, are highly polarized in excitation and are max- 
imally excited with electric fields that are oriented normal to the 
sample plane. At the highest surface coverage, the resonances of 
single hot spots are surprisingly narrow with widths as small as 
13 meV. The excitation spectra exhibit multiple, discrete resonances 
between 1.15-1.83 eV (1080-680 nm). Emission spectroscopy and 
super- resolution microscopy of the hot spots at their individual 
excitation resonances indicate that the emission at each resonance 
of one hot spot is localized to the same spatial region of the film 
within less than 10 nm. The supercontinuum generation appears to 
originate from a structural anomaly that is not readily identified with 
electron microscopy. In addition to providing a nearly background- 
free, non-invasive and non-perturbative probe of the plasmon reso- 
nances in these complex systems and unravelling their underlying 
structure-property relationships, nonlinear optical microscopy of 
single hot spots highlights the presence of extraordinarily high qual- 
ity factor, multi-resonant plasmonic excitations as was recently pre- 
dicted in evolutionary simulations of optical nanoantennae 18 . 

Methods 

Film growth. The growth process described by Saito et al. 27 was followed for the 
fabrication of Tollens films used in this work. Briefly, separate solutions of a silver- 
nitrate-based reagent and glucose were prepared in ultra-pure water. By combining 
the two solutions under ambient conditions (i.e., at room temperature and 
atmospheric pressure), the glucose reduces the reagent resulting in the deposition of 
silver on any glass surface in contact with the progressing reaction. Thus, by 
immersing a cleaned glass microscope coverslip into the solution, a silver film is 
deposited on the glass to form a Tollens substrate. The duration that the coverslip 
spends immersed in the reaction is defined as the growth time. 

Nonlinear excitation measurements. The excitation source (a Coherent Chameleon 
Ultra II Ti- Sapphire laser) produced excitation pulses with an average temporal width 
of 140 fs at an 80 MHz repetition rate and was tuneable between 1.82 eV and 
1.15 eV. The average spectral width of the excitation laser over the excitation range 
was 10 meV. The angle of incidence was 60°, and the orientation of the linear 
polarization was controlled with an achromatic XI 2 wave plate. Visible optical 
emission from the silver films was collected with a long working- distance objective 
(Olympus, 40X), passed through a short-pass emission filter that blocked light at 
energies below 1.90 eV, and imaged on a cooled charge-coupled device camera 
(Photometries Coolsnap HQ2) that was attached to an imaging spectrometer (PI 
Acton). All of the measurements were conducted in vacuum at room temperature. See 
the Supplementary Information for a list of the excitation densities used throughout 
this work. The FWHM of each resonance peak in the spectrum was determined by 
estimating the width of the peak at its half maximum using linear interpolation 
between neighbouring data points. The polarization modulation curves were 
analyzed by fitting a cos 4 curve to extract the values described in the text. 

Determination of in-plane dipole orientation. See the Supplementary Information. 

Filling factor calculation. See the Supplementary Information. 

Localization analysis. See the Supplementary Information. 
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